Saturated Ferromagnetism from Statistical Transmutation in Two Dimensions 
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The total spin of the ground state is calculated in the U — > do Hubbard model with uniform 
magnetic flux perpendicular to a square lattice, in the absence of Zeeman coupling. It is found that 
the saturated ferromagnetism emerges in a rather wide region in the space of the flux density <f> and 
the electron density n c . In particular, the saturated ferromagnetism at <f> = n e is induced by the 
formation of a spin- 1/2 boson, which is a composite of an electron and the unit flux quantum. 
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Ferromagnetism remains a challenging problem in spite 
of being among the best known phenomena in condensed 
matter physics. In particular, the origin of ferromag- 
netism in electron systems is fundamentally quantum- 
mechanical and non-perturbative Q. There are rather 
few established mechanisms of ferromagnetism, espe- 
cially those of the saturated (complete) ferromagnetism. 
Saturated ferromagnetism is defined as the ground state 
of the many-electron system with spin-independent inter- 
action having the maximum possible total spin. 

Nagaoka's theorem is one of few rigorous results on 
the saturated ferromagnetism |2| ■ It guarantees that the 
saturated-ferromagnetic state is the unique ground state 
when a single hole is inserted in the half-filled Hubbard 
model with infinite on-site repulsion U. Unfortunately, 
this theorem is limited to the single-hole case. Numerical 
studies suggest that Nagaoka ferromagnetism is unstable 
in the thermodynamic limit at finite hole densities 

The fiat-band ferromagnetism is another rigorous re- 
sult 0, Namely, the saturated ferromagnetism is 
proved rigorously under certain conditions, in electron 
systems with a (nearly) flat dispersion in the lowest band 
for a single electron. The ferromagnetism in a system 
with the low electron density and singular density of 
states near the Fermi level Q may also be related to the 
flat-band mechanism. However, as these results still have 
limited applicability, it is worth pursuing other mecha- 
nisms of (saturated) ferromagnetism. 

The difficulty in realizing the ferromagnetism in many- 
electron systems may be attributed to the Pauli principle 
for electrons. In the absence of interaction, the ground 
state of the system is generally paramagnetic, because 
the lower energy bands are filled with up and down spins. 
If we consider a system of bosons rather than fermions, 
the intrinsic tendency to favor paramagnetism may be 
absent. In fact, it was proved that in a continuous sys- 
tem with spinful bosons, one of the ground states is al- 
ways fully polarized if explicit spin-dependent interac- 
tions are absent [J Q ■ This statement holds also in a 
lattice model HEj- In particular, for the infinite-?/ Hub- 
bard model with spin- 1/2 bosons, the total spin of the 
ground state is shown to be maximal for all hole densi- 



ties, unlike in the electronic Hubbard model |9(. 

Thus, the saturated ferromagnetism could emerge in 
an electron system if the statistics of the electron is trans- 
muted to bosonic. In fact, the statistical transmutation is 
indeed possible in two-dimensional (2D) systems (lfUllj . 
and it has been applied to fractional quantum Hall ef- 
fect [13. 

Combining these ideas, we can expect that the 2D 
electron system in the presence of an external gauge 
(magnetic) field exhibits the (saturated) ferromagnetism 
thanks to the formation of the composite boson. Namely, 
when the applied magnetic field amounts to unit flux 
quantum per electron, the magnetic-flux quantum may 
be assigned to an electron. The composite particle con- 
sisting of an electron and the attached flux is then ex- 
pected to have spin 1/2 and to obey the Bose statistics 
with hard-core constraint. In this way, in the mean- field 
level, the original system can be mapped into a spin- 
1/2 boson system without magnetic field, which exhibits 
the saturated ferromagnetism 0,0, EJ However, as the 
"flux attachment" argument is not rigorous, whether this 
mechanism actually leads to the ferromagnetism in an 
electron system has to be checked. 

In quantum Hall systems in the continuum, the fully 
spin-polarized ground state is favored for the filling factor 
v = 1 (and in general, for v = \/rn with m odd), with- 
out the Zeeman energy [l3L HiL [HI . This is referred to 
as quantum Hall fcrromagnets |16| . While this ferromag- 
netism is usually associated with the antisymmetric na- 
ture of the orbital part of the wavefunction, it could also 
be regarded as a consequence of the formation of a spin- 
1/2 boson which is composed of an electron and m flux 
quanta |l6j . On the other hand, the dispersion in quan- 
tum Hall systems in the continuous space is completely 
flat (Landau levels). Thus the saturated ferromagnetism 
may also be understood as a special case of the flat-band 
ferromagnetism [J^- It is not clear whether the statis- 
tical transmutation is essential to realize the saturated 
ferromagnetism in quantum Hall systems. 

In order to clarify this issue, in this Letter we study an 
electron system on a lattice. We demonstrate an example 
of saturated ferromagnetism which is due entirely to sta- 
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tistical transmutation and is distinct from the fiat-band 
variety. 

Let us introduce the U — > oo Hubbard model on a 
square lattice, with the gauge (magnetic) flux per pla- 
quette [3]: 
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where (ij) refers to the nearest-neighbor pairs. Periodic 
boundary conditions are imposed in both directions, un- 
less explicitly mentioned otherwise. In the U — > oo limit, 
we have 
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where Cj CT = Ci a (l — m t -a\ which means that double 
occupancy at each site is excluded. We stress that, in 
our model J2J), we have not included the Zeeman coupling 
of spins to the magnetic field. The model is therefore 
completely isotropic in the spin space, and the total spin 
is a conserved quantum number. 

Exact numerical diagonalization for 4x4, -\/l8 x -\/l8, 
and \/20 x ^/20 clusters is employed in our study. We 
study the system with various values of the electron den- 
sity n e and the flux per plaquette 0. In particular, we 
need to investigate the case = n c where the statisti- 
cal transmutation to boson would occur. Under periodic 
boundary conditions, the total flux of the system is quan- 
tized to an integer. Thus can take only integral multi- 
ples of 1/TV, where TV is the number of sites (plaquettes). 
In order to study all the possible values of 0, we need 
to use the string gauge .19]: Choosing a plaquette S as 
a starting one, we draw TV — 1 outgoing arrows (strings) 
from the plaquette S, so that each plaquette other than S 
is the endpoint of a string. Then we set 0.y on a link ij to 
4>Mij taking account of the orientation, where Mij is the 
number of strings cutting the link ij. We have checked 
that the single-electron spectrum for a small system size 
(TV ~ 20) in the string gauge approximately reproduces 
the Hofstadter butterfly in the thermodynamic limit [2(| • 

The total spin Stat at zero temperature can be eval- 
uated from the expectation value of [Stat) 2 = (J2i Si) 2 
in the ground state, where Si is the spin operator at site 
I. In Fig. we show the scaled total spin Stot/Smax in 
the 0-n e plane. Here S max is given by TV e /2 with TV e be- 
ing the number of electrons. Red regions correspond to 
saturated-ferromagnetic states. In addition to Nagaoka 
ferromagnetism in the single- hole case with = 0, wo find 
two common features irrespective of the system size, (i) 
Saturated ferromagnetism appears along a straight line 
with = n c (or = 1 — n c ) except 0.6 n e 0.7. (ii) 



Saturated ferromagnetism appears in the triangular re- 
gion surrounded by three straight lines: n c = 0, = n e , 
and = 1 — n e . 

The result (i) confirms the expectation based on the 
statistical transmutation. Moreover, we find that, in 
most cases of = n c , the saturated ferromagnetism is ro- 
bust against twisting the boundary condition. This is in 
contrast to Nagaoka ferromagnetism, where the total spin 
of the ground state is c hang ed from maximum to zero as 
the boundary is twisted [21| . Besides, we have checked for 
TV = 16 and 18 that the saturated-ferromagnetic ground 
state at = n c is nondegenerate except for the trivial 
25'max + 1-fold degeneracy. This is consistent with the 
ferromagnetism for spin- 1/2 bosons 

In order to distinguish the ferromagnetism due to the 
statistical transmutation from possible "flat-band" vari- 
eties, we define the spectral functions 

= ^EK*"^ - l,TV i ;0)| C£T |vI/ o (TV T ,TV i ;0))| 5 
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x6[oj + E n (TVj - 1, TV i; 0) - £ (TV T , TV X ; 0) + fj], (3) 
= ^EK^W + l^i;0)|cl T |*o(TV T ,TV i; 0))| 2 

x5[w-E n (N r + l,N l] <l>) + E (N r ,N L ;^ + f i}. (4) 

Here /j, is the chemical potential, and (^(TVj, iVj.; 0)) 
denotes an eigenstate with energy E n (N^, TVj.; 0) in the 
system with TV| up-spins, TVj down-spins, and the flux 
0. We define the index n so that n = corresponds 
to the ground state with the given TVj and TVj. In the 
following, we set TV| = TVj_ = TV c /2 for TV G even. D ± 
can be estimated numerically by the continued-fraction 
method Below, we show the results for two values 
of n c ; n c = 4/20 and n c = 18/20 as representatives of 
the "low electron density" and "high electron density" 
regimes, respectively. 

First let us focus on the "low electron density" case, 
n c = 4/20. Figure |2Ia) shows the evolution of D ± (w) 
with varying 0. The sum of D + and D~ corresponds 
to the density of states (either occupied or unoccupied). 
Apparently it is always spread over a similar range of 
energy, representing the "bandwidth" which is about 8t 
although there is some 0-dependence. 

On the other hand, at this density, we find a crucial 
difference in the spectral function D~ corresponding to 
magnetism. Namely, D~ is concentrated in a narrow 
range of energy when the system exhibits the saturated 
ferromagnetism for = 4/20, 5/20, 9/20 (and 1 - 
0). In contrast, when the saturated ferromagnetism is 
absent, D~ is spread over a region of energy. Intuitively, 
the spectral function D~ corresponds to the density of 
states occupied by electrons. The narrow distribution of 
D~ compared to the "bandwidth" indicates a variant of 
the narrow or nearly flat-band ferromagnetism Q . 

In fact, there is more difference in D~ between the 
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(a)7V=16 (b)N=18 (c)N=20 




FIG. 1: Total spin of the ground state as functions of the electron density (n e ) and the flux per plaquette (<f>) in the 2D U — ► oo 
Hubbard model. Solid lines are straight lines with (f> = n c and (f> = 1 — n c , where statistical transmutation is expected. 



cases with and without saturated ferromagnetism, than 
what is visible in Fig. [21 D~ vanishes completely be- 
low a certain threshold (uj/t = -0.550,-0.229,-0.474 
respectively for <j> = 4/20,6/20 and 9/20) when the sys- 
tem exhibits the saturated ferromagnetism. In contrast, 
there is a continuous "shoulder" of low intensity (invis- 
ible in Fig. |2J down to much lower energy uj/t ~ —20 
when saturated ferromagnetism is absent. This seems to 
be consistent again with our interpretation. 

In particular, for <j> = 4/20(= n c ), the spectral func- 
tion D~ is localized within a narrow energy band below 
an apparent gap around the Fermi level. This appears 
similar to the quantum Hall ferromagnet at v = 1 in 
the continuum [16j. In this case, the statistical transmu- 
tation mechanism and the flat-band mechanism appear 
indistinguishable. 

Now let us discuss the saturated ferromagnetism ob- 
served in the "high electron density" regime, at <p = 
n c = 18/20. Figure |2Jb) shows the spectral functions 
in this case. Clearly, the spectral function D~ spreads 
over almost the entire bandwidth even though the sys- 
tem does exhibit the saturated ferromagnetism. Thus 
the ferromagnetism at <fi = n c = 18/20 is difficult to be 
understood in terms of the flat-band mechanism, and ap- 
pears to be exclusively due to the statistical transmuta- 
tion mechanism. In fact, changing the value of <j) destroys 
the saturated ferromagnetism at this electron density, as 
expected from the statistical transmutation scenario. 

In order to further confirm the statistical transmuta- 
tion scenario at <t> = n c , we define the following oper- 
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where Ji = 
In general, m 

denotes the number of magnetic-flux quanta, and we set 
m = 1 in the present case. On is the argument of the 
vector drawn from site i to site £. Note that the rela- 
tion 9u — Qu = ±1" holds for i ^ I. Then we can prove 
that these operators satisfy the commutation relations of 
bosons for two different sites. There is a large freedom 
in determining explicit values of du- Here we follow a 
prescription introduced in Ref. |23| . although this pre- 
scription unavoidably breaks translation invariance of a 
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FIG. 2: Spectral functions D ± in the 2D U -> oo Hubbard 
model, (a) n c = 4/20 with various values of <ft; (b) n e = 
18/20 with (j> = 18/20. The delta functions (vertical bars) 
are broadened by a Lorentzian with a width of O.li. The 
chemical potential is located at the zero energy. Cases with [F] 
exhibit saturated ferromagnetism. Colored regions represent 
_D~, while the dashed curves show D + . 



periodic cluster. The order parameter for the condensa- 
tion of the composite bosons may be defined by 

OB = ^EK*o( 7V T-l^i;^i)I^Tl*o(iV T ,^;0))| 2 . 



(5) 

We again set iVf = N± = N e /2 for 7V C even, and 
7V T = 7Vj_ + 1 = (N c + 1) /2 for N c odd. FigureEta) shows 
Ob as a function of n e (= 0). The order parameter has a 
pronounced enhancement in both the high-density region 
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FIG. 3: Order parameter of spin- 1/2 bosons in the 2D U — » oo 
Hubbard model, (a) n c (= <ji)-dependence; (b) ^-dependence 
for n c = 6/20; (c) (^-dependence for n c = 16/20. In (b) and 
(c), data are not plotted where the initial and/or final states 
in Eq. © have degeneracy. 



to observe the magnetic-field effect at a modest magnetic 
field of a few tesla. In this situation, the orbital motion 
rather than the Zeeman effect could be essential for the 
emergence of ferromagnctism. Another possibility is to 
induce the effective gau ge field internally without an ap- 
plied magnetic field [26j. 
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(0.7 < n c < 1) and the low-density one (0 < n c 0.3). 
Figures Efb) andEJc) depict the (^-dependence of Ob for 
n c = 6/20 and 16/20, respectively. We find a salient 
growth at = n e for both densities. We have also ob- 
served for n c = 6/20 that the order parameter given by 
Eq. JSJ with m = — 1 has a peak at = 1 — n . These 
results again support the ferromagnetism based on the 
statistical transmutation. 

Finally, we discuss the region 0.6 n c 0.7, where 
the saturated ferromagnetism is absent in spite of <j> = n e . 
As a candidate of the competing order, we consider 
the spin chirality defined by the order parameter 
Xch = (l/N)Y,i(Si ■ S e+ y x Se+x), where x (y) is the 
unit vector along x (y) direction, and (• ■ ■) denotes the 
expectation value in the ground state. In fact, Nagaoka 
ferromagnetism in the single-hole case is known to be de- 
stroyed by development of the spin chir ality in the pres- 
ence of a perpendicular magnetic field [lq. Along the 
line <j> = n e , we have confirmed that the spin chirality 
vanishes when the system exhibits saturated ferromag- 
netism. On the other hand, the chiral order is indeed 
developed when the saturated ferromagnetism is absent 
(not shown). 

In summary, we have calculated the total spin of the 
ground state in the U — > oo Hubbard model with mag- 
netic flux (0) perpendicular to a square lattice and re- 
vealed regions of saturated ferromagnctism. The satu- 
rated ferromagnetism at cf> = n c is argued to be due to 
formation of spinful composite bosons. Statistical trans- 
mutation may therefore play a key role in ferromagnetism 
in strongly correlated systems, just as it did in fractional 
quantum Hall effect. 

The present mechanism may be relevant to future ex- 
periments on an artificial crystal of a square lattice with 
quantum dots (i.e., a quantum dot superlattice) |2fij . A 
large lattice constant of such a crystal would enable us 
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